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Abstract of correspondent: US61 75725 
A receiver of modulated carrier-wave signals 
equipped with a pseudo-random code time 
marker, particularly for radio navigation, 
comprises coincidence modules each of which 
comprises chronometric memory medja which 
receives a wanted signal, provides storage as 
data of the time offset of the code offset 
estimates and the frequency and phase carrier- 
wave status, and on that basis provides, with a 
local clock, a local carrier-wave image and that 
least one local code repeat for the relevant signal 
and code and carrier-wave slaving control 
functions providing correlation of the wanted 
signal and the local code repeat on the basis of a 
frequency and phase offset signal between the 
carrier-wave of the wanted signal and the carrier- 
wave image. It further comprises Fourier 
transformation media for reception of the offset 
signal and management/decision-making media 
which provide invalidation of time offset data 
where the Fourier transformation indicates the 
presence of energy outside of the vicinity of a 
central coincidence line 
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(54) Recepteur de radionavigation, en particulier du type GPS 



(57) Un recepteur de signaux a porteuse modulee 
munie d'un marqueur temporel a code pseudo-a!6atoi- 
re, notamment pour la radionavigation, comprend des 
modules a coincidence comportant chacun des moyens 
de memoire chronometrique (71-75,60-89) qui recep- 
tionne un signal utile, stockent comme donnees de ca- 
lage temporel des estimees de calage de code (C_PHI) 
et d'etat de porteuse en frequence (P_DF) et phase 
(P_PHI), et en tirent, avec une horloge locale, une ima- 
ge locale de porteuse et au moins une replique locale 



du code, pour le signal concerne, et des asservisse- 
ments de code et de porteuse (40,50,59), correlant le 
signal utile et la replique locale de code a i'aide d'un 
signal d'§cart de frequence et de phase (DVCa) entre 
la porteuse du signal utile et Pimage de porteuse. II com- 
prend en outre des moyens de transformation de Fou- 
rier recevant le signal d'6cart (DVCa) et des moyens de 
gestion/decision (90) qui invalident les donnees de ca- 
lage temporel lorsque la transformation de Fourier indi- 
que la presence d'energie en dehors du voisinage d'une 
raie centrale de coincidence. 
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D scription 

[0001] L'invention concerne la radionavigation, qui 
utilise le temps de propagation d'une onde electroma- 
gnetique & marqueur temporel entre un emetteur et un 
recepteur. 

[0002] Dans les systemes de navigation recents, le 
marqueur temporel de I'onde porteuse est du type pseu- 
do-bruit repetitif; on utilise en pratique un code pseudo- 
aleatoire. A la reception, le temps de propagation de 
I'onde se manifeste a la fois par un decalage temporel 
du code pseudo-aleatoire et par un decalage de phase 
de la porteuse. En presence d'un mouvement relatif sur 
I'axe emetteur/recepteur (ligne de visee), il s'y ajoute un 
decalage de frequence de la porteuse, du k I'effet Dop- 
pler. 

[0003] Differentes sources d'erreur existent: certai- 
nes, liees a la traversee de la troposphere ou de I'ionos- 
phere, ou bien & la variation de conductivity des sols, 
affectent de maniere generate la vitesse de propagation 
des ondes; d'autres sources d'erreurs tiennent au fait 
que le trajet parcouru par I'onde radioelectrique jusqu'au 
recepteur n'est pas rectiligne, d'autres encore au fait 
que le recepteur pergoit une combinaison de ditfe>ents 
trajets issus de reflexions diverses, generalement dues 
a des surfaces voisines de I'antenne de reception. C'est 
ce qu'on appelle les erreurs de trajets multiples: au trajet 
direct (le plus court) se superposent d'autres trajets in- 
desires. 

[0004] Dans les systemes de radionavigation les plus 
utilises a ce jour, dits GPS et GLONASS, les emissions 
proviennent de satellites. 

[0005] Actuellement, le probleme principal est de 
combattre les effets des trajets multiples, ainsi qu'en t6- 
moigne I'article "Conquering multipath: The GPS accu- 
racy battle", Lawrence R. WEILL, revue GPS World, 
Avril 1997, pp 59-66. Tel est le but principal de la pr6- 
sente invention. 

[0006] L'invention s'applique & un dispositif recepteur 
de signaux ou ondes radio-6lectriques, qui portent un 
marqueur temporel a pseudo-bruit repetitif, du genre co- 
de pseudo-aleatoire : en particulier pour la radionaviga- 
tion. 

[0007] Un tel dispositif comporte une reception haute 
frequence, dont la sortie (ou "signal utile") est distribute 
a plusieurs canaux ou modules, dont chacun estaffecte 
& la recherche d'une coincidence (calage temporel) 
av c un signal radio-electrique choisi respectif, distin- 
gue par son code et/ou par sa frequence porteuse. II s'y 
ajoute des moyens d'exploitation des donnees de cala- 
ge temporel fournies par ces modules a recherche de 
coincidence (ou plus simplement : module & coinciden- 
ce). 

[0008] Chaque module & coincidence comprend 
d'une part des moyens forrhant memoire chronometri- 
qu , aptes a recevoir et/ou stocker une estimee de ca- 
lage de code et une estimee d'etat de porteuse en fre- 
quence et phase, ainsi qu'a en tirer, a I'aide d'une hor- 



loge locale, une image locale de porteuse et au moins 
une replique locale du code, pour le signal radio-elec- 
trique concerne. II comprend d'autre part des asservis- 
sements visant a faire coincider la replique locale du co- 

s de avec le code recu. Ces asservissements compren- 
nent un asservissement de code, operant par correla- 
tion entre le signal utile et la replique locale de code, et 
un asservissement de porteuse, qui travaille & I'aide 
d'un signal d'ecart de frequence et de phase entre la 

10 porteuse presente dans le signal utile et ladite image de 
porteuse. 

[0009] Un tel module a coincidence peut fournir une 
indication sur le bon fonctionnement de ses asservisse- 
ments. On sait en tenir compte d'une part pour le choix 
is de signaux radio-electriques incidents qui donnent des 
resultats satisfaisants, d'autre part au niveau de I'exploi- 
tation des donnees de calage temporel. 
[0010] Selon un aspect de l'invention, le dispositif re- 
cepteur comprend en outre des moyens capables d'et- 
20 fecluer une transformation de Fourier recevant un signal 
repr§sentatif de I'ecart de frequence et phase de por- 
teuse. A partir de la transformed de Fourier, des moyens 
de decision peuvent invaliderau moins partiellement les 
donnees de calage temporel, en particulier lorsque la 
2S transformation de Fourier indique la presence d'energie 
en dehors du voisinage d'une raie centrale, cette der- 
niere correspondant k la coincidence. 
[0011] On precisera plus loin differentes caractdristi- 
ques plus delaillees du fonctionnement de cette sur- 
30 vetllance frequentielle. 

[0012] Dans un mode de realisation pr6f6rentie! du 
module a coincidence, les moyens formant memoire 
chronometrique comprennent un oscillateur de porteu- 
se pour definir I'image locale de porteuse, un oscillateur 
35 de code, et un generateur d'au moins une replique lo- 
cale du code. De leur cote ; les asservissements de code 
et de porteuse comprennent une voie de correlation 
ponctuelle, comportant une demodulation selon ladite 
image de porteuse, et fournissant un signal de coinci- 
de dence entre la position temporelle du code dans le si- 
gnal utile et son estimee, et au moins une voie de cor- 
relation differentielle, comportant une demodulation se- 
lon ladite image de porteuse, et fournissant un signal 
d'ecart entre la position temporelle du code dans le si- 
45 gnal utile et son estimee. Le signal de coincidence sert 
de grandeur d'ecart pour une boucle de porteuse a ver- 
rouillage de phase, laquelle pilote le g6nerateurd'image 
locale de porteuse; et le signal d'ecart de code sert de 
grandeur d'ecart pour une boucle de code a verrouillage 
so de retard, dont la sortie est combinee a celle de la boucle 
de porteuse a verrouillage de phase, pour piloter le ge- 
nerateur de repliques de code. 

[0013] Un autre aspect de l'invention, tres interessant 
en lui-nrTme, est I'analys© "temporelle". Selon celle-ci, 
ss on etablit une autre indication de la presence de multi- 
ples trajets en utilisant le signal d'ecart issu d'au moins 
une voie de correlation differentielle. 
[0014] Pour mettre en oeuvre I'analyse temporelle, il 
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est actuellement prolere que chaque module a coinci- 
dence comprenne an moms une seconde vole de cor- 
relation differentiellc. possedant un espacement diffe- 
rent de celui de la premiere On prevoit alors des 
moyens pour utiliser le signal d'ecart issu de cette se- 
conde voie de correlation, seul et/ou par sa compari- 
son (difference) avec le premier signal d'ecart, en tant 
qu'indication de la presence de multiples trajets. 
[0015] On precisera plus loin differentes caracteristi- 
ques plus detaillees du fonctionnement de cette sur- 
veillance ou analyse temporelle. 
[0016] D'autres caractenstiques et avantages de in- 
vention apparaitront a I'examen de la description de- 
taillee ci-apres, et des dessms annexes, sur lesquels: 

la figure 1 est le schema do pnncipe general d'un 
recepteur GPS; 

la figure 2 est le schema plus delailte d'un mode de 
realisation particulier de recepteur GPS, auquel in- 
vention peut s'appliquer. 

les figures 3 et 4 illusti ent deux details du recepteur 
de la figure 2; 

la figure 5 est un schema equivalent permettant de 
mieux comprendre la generation des repliques lo- 
cales appliquees a des corrclatcurs dans le recep- 
teur des figures 2 a 4; 

la figure 6 est une vue partielle illustrant Tune des 
caracteristiques de la presente invention; 
la figure 7 illustre le mecanisme associe a la sur- 
veillance Irequentielle selon invention; 
les figures 8A-1 a 8G-1 sont des diagrammes fre- 
quentiels illustrant differentes situations de multi- 
ples trajets, speculates et/ou diff us; 
les figures 8A-2 a 8G-2 sont des diagrammes de 
transformees de Fourier, equivalents respective- 
ment aux diagrammes des figures 8A-1 a 8G-1; 
les figures 9A a 9C sont trois diagrammes tempo- 
rels illustrant respectivement les sorties des deux 
voies de correlation differentielle B et C du recep- 
teur des figures 2 a 4, ayant des espacements dif- 
ferents, et le signal difference entre ces deux 
sorties ; et 

la figure 10 illustre schematiquement les moyens 
de base de I'analyse temporelle selon I' invention. 

[0017] Les dessins annexes sont, pour I'essenttel, de 
caractere certain. En consequence, ils pourront non 
seulement servir a mieux faire comprendre la descrip- 
tion, mais aussi contribuer a la definition de invention, 
le cas echeant. 

[0018] La radionavigation utilise surtout les codes 
pseudo-aleatoires, bien que Ton puisse concevoir 
d'autr s types de marqueurs temporels assimilables a 
un bruit. Un code pseudo-aleatoire st un enchatne- 
ment repetitif d'une sequence pseudo-aleatoire. On 
f'appelle aussi code a bruit pseudo-aleatoire, ou bien co- 
de PRN (de I'anglais "Pseudo Random Noise"). 
[0019] Une facon d'engendrer une sequence pseudo- 



aleatoire est decrite dans FR-A-2 248 517, notamment 
en reference a ses figures 1 a 3. Une sequence pseudo- 
aleatoire peut etre consideree comme une suite de bits, 
parcourue avec une periode d'horloge bien definie. L'ex- 

s pression appelee "moment de code", ou plus frequem- 
ment "chip", designe un bit de la sequence, et, par ex- 
tension, la duree (Dmc) d'un tel bit. 
[0020] Une sequence pseudo-aleatoire possede la 
propriete suivante: sa fonction d'autocor relation est nul- 

70 le, sauf au voisinage du decalage temporel nul, ou elle 
prend une allure triangulaire, avec un maximum en cas 
de synchronisme. La largeur de la base du triangle est 
de +/- Dmc (voir la figure 13A de FR-A-2 248 517 et la 
description correspondante). La valeur zero indique un 

15 niveau de rejection qui depend de la longueur du code. 
II en decoule que Ton peut engendrer une variete de se- 
quences pseudo-aleatoires de meme duree ou "epo- 
que" (Dprn), et dont les fonction s d'intercor relation sont 
nulles. 

20 [0021] Dans le systeme GPS, differents satellites 
emettent des porteuses, a une frequence d'environ 1 
gigahertz, modulees par des codes pseudo-aleatoires 
respectifs, et specifiques a la fois du satellite et de la 
frequence porteuse. Plus precisement, un satellite don- 

2B no emet actuellement deux porteuses, munies de codes 
PRN respectifs notes 'C/A" (pour "Coarse Acquisition" 
ou "Clear Acquisition"), et "P" (pour "Precise"). Le pre- 
mier code C/A est public, alors que le second code note 
P, presente une version notee P(Y), intrinsequement 

30 plus precise que le code C/A, mais soumise a un en- 
cryptage qui reste confidentiel. La polarisation de I'onde 
porteuse est circulate droite. 

[0022] Comme les recepteurs GPS d'aujourd'hui uti- 
iisent surtout le code "C/A", non brouille, la suite de cette 
35 description est limitee a ce code et a la porteuse corres- 
pondante, a 1575,42 MHz, qui comporte deux 
modulations : 

une modulation numerique rapide par inversion de 
40 phase (BPSK pour Bi-Phase Shift Keying) selon un 
code pseudo-aleatoire specifique de 1023 bits, pul- 
se a 1,023 MHz, avec Dprn = 1 ms, et 
une modulation numerique beaucoup plus lente (50 
Hz) definissant des donnees d'etat du satellite, que 
« |'on appelle "message de navigation". 

[0023] II est clair que invention ne se iimite pas a rap- 
plication au seul code C/A. 

[0024] Les satellites GPS ont une reference de temps 
so commune. Le message de navigation d'un satellite 
("ephemerides" et "almanachs") comprend notamment 
la derive de son horloge propre (par rapport a une refe- 
rence de temps), et les elements necessaires pour cal- 
culer sa position (parametr s d'orbite). En bref , les pa- 
ss rametres de leur orbite etant connus, les satellites sont 
a considerer comme des postes emett urs dont les po- 
sitions sont determinates. 

[0025] Le schema synoptique general d'un recepteur 
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GPS est donne sur la figure 1. A partir des ondes elec- 
tromagnet iques recues de differents satellites GPS (si- 
gnaux bruts), le bloc 1001 realise la reception propre- 
ment dite, qui convertit les ondes en signaux electri- 
ques. Viennent ensuite un bloc fonctionnel 1002 de con- 
version de frequence, puis un bloc fonctionnel 1 003, qui 
correspond a I'echantiUonnage et a la quantification. 
[0026] Le bloc 1 004 correspond aux correlations, ain- 
si qu'aux autres mesures de temps et demodulations 
effectuees sur la porteuse. A ce stade. une correlation 
du code recu avec la replique locale du generateur de 
code correspondant au satellite a acquerir (engendree 
sur la base d'une horloge locale) va permettre de deter- 
miner le temps d'arrivee de I'onde recue, et de retrouver 
la modulation lente qu'elle porte. En fait, si les deux co- 
des d'entree du correlateur sont en coincidence ou 
presque (a I'interieur de la reponse triangulaire preci- 
tee), on obtient une tension en sortie; sinon, la sortie du 
correlateur est nulle, c'est a dire du bruit. Et c'est en ac- 
cumulant cetle tension pendant un certain lemps que 
Ton obtient une information significative sur la coinci- 
dence (compte-tenu du fait que le signal recu est au des- 
sous du bruit dans la bande d'entree du recepteur), et. 
que Ton peut retrouver le message de navigation. Com- 
ma Poperationde correlation se trouve ainsi otalee dans 
le temps, de tels systemes se contentent d'une puissan- 
ce d'emission faible, correspondant a un niveau de re- 
ception situe au dessous du niveau de bruit. 
[0027] En pratique, le bloc 1004 comporte deux bou- 
cles d'asservissement. Une boucle d'asservissement 
de porteuse "poursuit" la frequence et la phase de la 
porteuse, telle que regue. Une boucle d'asservissement 
de code "poursuit 1 ' le synchronisme entre le code recu 
et la replique locale correspondante. 
[0028] La boucle de porteuse memorise done la pha- 
se totale estimee <X>j (phi majuscule) de la porteuse. L'ex- 
pression "phase totale" comprend ici le decalage de 
phase subi par la porteuse pendant le temps de trajet 
des ondes et I'eff et Doppler du au mouvement relatif sa- 
tellite/recepteur. En fait, il est plus commode de separer 
la phase totale en un ecart de frequence f it et une phase 
(|>j (phi minuscule), conformement au fonctionnement 
d'une boucle de porteuse du second ordre, qui realise 
d'abord I'accrochage en frequence, puis I'accrochage 
en phase. L'ecart de frequence fj correspond a un de- 
calage Doppler, done a une Vitesse relative; c'est pour- 
quoi on I'appelle aussi "Vitesse". 

[0029] La boucle de code memorise un retard estime 
du code t j, par rapport a Phorloge locale. L'ambiguite\ 
egale a une longueur de code (Dprn = 1 ms, soit environ 
300 km), est levee a Paide du message de navigation, 
et eventuellement d'une estimee de position disponible 
par ailleurs. 

[0030] De facon connue, le recept ur peut alors pren- 
dre differents etats, ce qui correspond a I'unite fonction- 
nelle detest 11 00 dite "aiguillagesuivant I'etatdela ligne 
de visee". Dans une phase d'initialisation complete, on 
realisera d'abord en 1105 I'accrochage ou acquisition. 



Ceci fait, on pourra passer a la phase de poursuite 1 1 06. 
Par contre, si un accrochage a deja ete etabli, et que 
Ton est deja passe anterieurement en phase de pour- 
suite, mais que le synchronisme a ete perdu, on pourra 
5 afors realiser un reaccrochage comme indiqud en 11 07. 
Ce reaccrochage peut tenir compte de donnees en pro- 
venance de capteurs exterieurs au recepteur GPS, 
comme indique en 1110. 

[0031] Lorsqu'on est en phase de poursuite 1106, il 
10 est possible de realiser le decodage des messages de 
navigation (bloc 1108). En 1119 survient un test de co- 
herence de I'ensemble des mesures, qui peut tenir 
compte aussi de donnees en provenance des capteurs 
exterieurs au recepteur (1110). Ce test de coherence 
*5 fait avantageusement intervenir les techniques dites 
RAIM (Receiver Autonomous Integrity Monitoring). Ces 
techniques appliquees dans le cas de mesures redon- 
dantes, ont pourobjet de determiner les lignes de visee 
qui risquent d'entrainer une erreur de position, et de les 
20 ecarler du traitement de determination de position. Elles 
sont dec rites notamment dans le Journal of the Institute 
of Navigation : 

vol.35 N°4, hiver 88-89, article "Navigation System 
2S Integrity Monitoring Using Redundant Measure- 
ments" de Mark A. STURZA, 
vol.35 N°2, ete 88, article "Autonomous GPS Inte- 
grity Monitoring Using the Pseudorange Residual" 
de BRADFORD W. PARKINSON et Penina 
30 AXELRAD, et 

vol.39 N°3, automne 92, article "A Baseline GPS 
RAIM Scheme and a Note on the Equivalence of 
Three RAIM Methods" de R. GROVER BROWN. 

35 [0032] On en trouve egalement une description dans 
I'Appendix "O" de "Minimum Operational Performance 
of Standards for GPS/Wide Area Augmentation System 
Airborne Equipement" : RTCA, Doc N° RTCA/DO-229, 
January 16,1996, prepared by SC-159, dont I'Appendix 

40 *E " donne les references d'autres articles. 

[0033] A partir des instants de reception des ondes 
obtenus en 1001, le bloc 1120 applique des techniques 
du genre triangulation, qui permettent de determiner la 
position et la Vitesse du recepteur GPS. On obtient done 

45 un vecteur de position P et un vecteur Vitesse V. Les 
coordonnees de position sont definies par rapport a un 
ellipsoide de reference dit WGS 84 (World Geodetic 
System), et different du geoTde. 

[0034] Plusieurs mesures de temps de trajet (au 
50 moins quatre) sont necessaires pour determiner la po- 
sition du recepteur On travaille en principe avec redon- 
dance, pour augmenter la precision, et supporter plus 
facilement les consequences de la perte de contact 
avec run des satellit s. Pour determiner la position du 
55 recepteur, le bloc 1 1 20 doit egalement calculer une es- 
timation du "biais d'horloge", c'est a dire du decalage de 
temps entre I'horloge du recepteur et le "temps GPS", . 
a savoir le temps general (de reference) du systeme 
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commun a tous les satellites. 

[0035] On dispose deja de certains moyens pour ten- 
ter de s'affranchir des signaux parasites que creent les 
multiples trajets. 

[0036] Une premiere classe de moyens concerne les 
aeriens. Une bonne rejection des signaux parasites en 
polarisation inverse peut etre procuree par des anten- 
nes de type helicol'dal, qui attenuent la puissance des 
multiples trajets ayant subi un nombre impair de re- 
flexions. Des plans masses ou des plans absorbants 
peuvent etre disposes sous les antennes pour attenuer 
Peffet des reflexions venant du sol, ou bien dans une 
situation d'incidence rasante. On connait aussi des an- 
tennes a structure annulaire, plus precisement denom- 
mees "choke ring", et qui constituent des systemes al- 
veolaires formant pieges. On peut encore utiliser des 
antennes a formation de faisceaux, a diagramme directif 
orientable, afin d'ecarter des multiples trajets Indesires, 
ou du moins de les att6nuer. Sont encore connus des 
systemes multi-antennes qui se fondent sur la decorre- 
lation spatiale des multiples trajets, afin d'eliminer ceux- 
ci. 

[0037] Ces differentes solutions peuvent etre mises 
en oeuvre dans les dispositifs ici decrits, mais ne cons- 
tituent pas en eux-memes le coeur de la presente in- 
vention. 

[0038] On cherche egalement a reduire I'effet des 
multiples trajets au niveau de la conception des recep- 
tees eux-memes. 

[0039] Le procede dit "harmonique" exploite la diver- 
gence code/phase observee sur une ligne de visee (la 
ligne de visee est la ligne qui joint la ligne du trajet effectif 
le plus court entre un satellite et le recepteur). Cette 
technique, mise en oeuvre par la Demanderesse, est 
fondee sur le tait que I'erreur observee sur la porteuse 
est plus faible que celle que Ton observe sur le code. 
En procedant a un lissage du code par la phase, on peut 
observer I'ecart code-phase, et en deduire la presence 
des multiples trajets. En effet, une augmentation de la 
moyenne de I'ecart code-phase observee sur une ligne 
de visee est symptomatique de la presence de multiples 
trajets, ce qui peut d'ailleurs servir a titre comp!6men- 
taire dans la mise en oeuvre de la presente invention. 
[0040] Le document FR-A-2 698 966 decrit une autre 
technique consistant a analyser le front avant du pic de 
correlation, remarque etant faite que les multiples trajets 
(physiques) sont toujours en retard sur le trajet direct, 
et que le front avant est peu ou pas d6forme. On obtient 
ainsi une amelioration du comportement du recepteur. 
[0041] II a 6galement 6te propose d' exploiter la cour- 
be de correlation, par analyse de la variation de sa pen- 
te, et par utilisation de multiples correlateurs s'appuyant 
sur la deformation de cette courbe par rapport a sa for- 
me canonique. Cett technique, decrite notamment 
dans le Brevet US-A-5 390 207, donne des resultats in- 
teressants lorsqu'un nombre restreint de multiples tra- 
jets existe, avec un etat suffisamment stationnaire. 
[0042] Par ailleurs, il a 6galement ete propose de re- 



duire I'espacement de la correlation differentielle dans 
une boucle DLL qui travaille en phase de poursuite. 
comme indique par exemple dans le fascicule de Brevet 
EP-A-0 488 739. 
5 [0043] Les proc6des ci-dessus visent a detecter les 
trajets multiples, et, pour certains, a s'en affranchir par- 
tiellement. 

[0044] La Demanderesse a observe que, si ces pro- 
cedes connus conviennent avec certaines classes bien 

10 repertories de trajets multiples, par contre ils ne sont 
pas d'application generate, et risquent meme dans cer- 
tains cas d'avoir un effet oppose a celui recherche, en 
presence de trajets multiples qui n'appartiennent pas 
aux classes traitees avec succes. A cot6 de cela, en 

15 multipliant les correlateurs (qui sont a multiplier en outre 
par le nombre de satellites traites), on complique con- 
siderablement le recepteur. 

[0045] En consequence, apres avoir etudie en profon- 
deur le phenomena en cause, la Demanderesse a pre- 

20 fere se focaliser sur des amenagemenls du recepteur 
qui permettent la detection de trajets multiples, afin 
d'empecher le recepteur d'utiliser des mesures erro- 
nees pour se localiser (ou bien d'utiliser ces mesures 
erronees a d'autres fins). La redondance des mesures 

25 disponibles et la probability faible d'observer des trajets 
multiples simultanement sur plusieurs lignes de visee 
permettent done, avec un investissement materiel mini- 
me, d'accroTtre de fagon significative la quality du posi- 
tionnement. 

30 [0046] Bien que I'invention puisse s'appliquer a diffe- 
rentes structures de recepteurs GPS, elle s'harmonise 
particulierement bien avec un recepteur a trois voies de 
correlation, que Ton decrira des maintenant. La figure 2 
est un schema plus detaille d'un tel recepteur, limite a 

35 un canal de satellite pour simplifier. 

[0047] Le recepteur commence par une antenne 101 
et un 6tage de haute frequence 102 (bloc 1001 de la 
figure 1). Vient ensuite un amplificateur 103, 1'ensemble 
operant sur une bande qui est ici de 20 MHz, a une fre- 

40 quence de 1575,42 MHz, dans I'exemple choisi. De fa- 
gon generale, ia largeur de bande en haute frequence 
est comprise entre environ 2 Mhz et quelques dizaines 
de MHz. 

[0048] Un melangeur 105 regoit la sortie de I'amplifi- 
45 cateur 103, ainsi qu'un signal local de 1400 MHz d'un 
oscillateur local 107; il fournit en sortie un signal utile 
analogique d'une frequence de 175,42 MHz. Ce signal 
est applique a un filtre passe-bande 109, centre sur 
175,42 MHz, avec une r§ponse de +/- 10 MHz a -3 dB, 
50 et +/- 25 MHz a -15 dB. Le bloc 1002 de la figure 1 cor- 
respond a ces organes 104 a 109, auxquels on peut 
ajouter I'amplificateur 103. 

[0049] Le filtre 109 est suivi d'un amplificateur 110, 
puis d'un convertiss ur analogique num6rique 112, dont 
55 ia sortie st le signal utile, cette fois sous forme nume- 
rique. L'amplificateur 110 est de preference a gain va- 
riable, son gain etant commande en fonction du niveau 
de signal a la conversion, de fagon a optimiser le nom- 
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bre de chiffres signil-calils en sortie du convertisseur (2 
bits). 

[0050] Un generateur dc signaux d'horloge 20 produit 
les differents signaux dhoMoges utiles au recepteur, et 
qui sont, sauf mention conlraire tous lies entre eux. Ain- 
si ce generateur 20 fournit par exemple une frequence 
de 20 MHz, qui sert a piloter (par multiplication) Poscilla- 
teur local 107. Ce dernier comporte un oscillateur com- 
mande en tension (VCO) operant a 20 MHz, et muni 
d'une boucle a verrouillage de phase (PLL) sur cette fre- 
quence pilote de 20 MHz A partir de la sortie de I'os- 
cillateur, des generateurs de frequence fournissent la 
frequence de 100 MHz qui cadence Pechantillonnage 
au sein du convertisseur analogique numerique 112, et 
les autres frequences utiles de 1400 MHz et 25 MHz, 
essentiellement. 

[0051] Le signal analogique possede une frequence 
de 175,42 MHz. II est sous-echant»llonne a 100 MHz, ce 
qui equivaut a une transposition de frequence a 200 
MHz. La sortie du convertisseur analogique-numerique 
112 peut done etre consideiee comme un signal nume- 
rique centre sur 24,58 MHz. Exprime sur 2 bits, ce signal 
numerique, qui comporte ici 4 echantillons par periode, 
est applique a un echantillonneurcomplexe 120, detaille 
sur la figure 3. 

[0052] En tete, I'echantillonneur comporte en 121 un 
transposeur de frequence a bande laterale unique de 
420 Kilohertz, ramenant la frequence de travail a exac- 
tement 25 MHz (24,58+0,42). Cadence en sortie par la 
frequence H100M (le meme signal d'horloge que pour 
le convertisseur 112), le signal issu de ce bloc 121 est 
applique a deux operateurs 1 23I et 1 23Q cadences res- 
pectivement par des horloges H25M-I et H25M-Q, a 25 
MHz, en quadrature Tune de Pautre. Ces operateurs 1 23 
sont agences pour fournir des signaux numeriques en 
bande de base, exprimes sur 4 bits, a savoir respecti- 
vement un signal en phase (sortie I), et un signal en qua- 
drature (sortie Q). Par analogie avec les nombres com- 
plexes, on parle aussi des parties reelle I et imaginaire 
Q. 

[0053] Le bloc 1003 de la figure 1 comprend Pensem- 
ble qui va du convertisseur analogique/numerique 112, 
precede de son amplificateur a gain variable 110, jus- 
qu'a I'echantillonneur numerique 120. 
[0054] Les deux sorties I et Q de I'echantillonneur nu- 
merique complexe 120 sont distributes a trois votes 
(bloc 1004 de la figure 1). Ces trois voies sont distin- 
guees par les letlres A, B et C, en correspondance des 
trois correlations egalement notees A, B et C, auxquel- 
les sont associes trois decalages possibles du code 
pseudo-aleatoire determinant le satellite a poursuivre. 
On verra plus loin que la voie B est •differentielle", e'est 
a dire qu'elle travaille sur la difference entre deux repli- 
quesd'un meme code pseudo-aleatoire, espaceestem- 
porellement Tune de I'autre autour du decalage associe 
a cette voie B. La voie C peut etre normale (code C1), 
ou bien differ ntielle (code C2) 

[0055] La figure 4 detaille a titre d'exemple la voi A, 



qui apparait en haut sur la figure 2. Les signaux en pha- 
se I et en quadrature Q sont appliques respectrvement 
a deux correlateurs numeriques 31 1 A et 31 2A, recevant 
par ailleurs le code A, et dont les sorties vont conjointe- 

5 ment vers un demodutateur de porteuse 32A. 

[0056] Pour le signal I (partie reelle), la sortie du cor- 
relateur 311 A (4 bits) est appliquee d'abord a un som- 
mateur 321, qui fait un cumul sur 250 echantillons nu- 
meriques (9 bits). Ce cumul est applique a deux multi- 

10 plieurs numeriques 322 et 323, recevant par ailleurs res- 
-pectivement les informations COSinus et SINus de la 
phase de porteuse (P_PHI), exprimees chacune sur 6 
bits. Les resultats respectifs, notes Icos et Isin, sont ex- 
primes sur 9 bits. Le schema est le meme pour le signal 

is Q (partie imaginaire), avec le correlateur 31 2A (4 bits), 
le sommateur 331, et les deux multiplieurs numeriques 
332 et 333, dont les resultats respectifs, notes Qcos et 
Qsin, sont exprimes sur 9 bits. 

[0057] A partir de la, les sorties principals de la voie 
20 a, en phase et en quadrature, comprennent un signal 
LA, fourni par un additionneur numerique 325 qui cons- 
truit Icos + Qsin : sur 10 bits, ce qu'un sommateur 327 
accumule sur 16 bits, et un signal Q_A, fourni par un 
soustracteur numerique 335 qui construit Qcos - Isin, 
25 sur 10 bits, ce qu'un sommateur 337 accumule sur 16 
bits. 

[0058] Les accumulations en 327 et 337 se font sur 
100 ou 200 echantillons, de facon commutable, ce qui 
correspond a un temps d'integration de 1 ou 2 millise- 

30 condes (base de temps a 1 kHz ou 500 Hz). En effet, 
partant d'environ 25 MHz, on a realise des cumuls sur 
250 fois 100 ou 200 echantillons numeriques. 
[0059] Les deux autres voies de la figure 2 sont de 
structure identique a ce qui vient d'dtre decrit, et portent 

35 les memes references numeriques, mais avec les suf- 
fixes B et C, respectivement, la ou les correlateurs re- 
coivent le code B et le code C. 

[0060] En outre, pour la voie A seulement (trait tirete 
sur ia figure 4), la sortie de Padditionneur 325 est appli- 

40 quee a un sommateur 341 , operant sur 50 echantillons. 
Le bit de signe de la sortie du sommateur 341 , disponi- 
ble toutes les 0,5 ms, est memorise et serialise dans la 
memoire 342 pour fournir 2 ou 4 informations successi- 
ves de signe, selon que I'on travaille a 1 kHz (1 ms) ou 

45 500 Hz (2 ms). Cette information de signe, notee 
SIGNJA_2K, correspond a la modulation lente porteu- 
se des donnees d'etat, "ephemeride" et "almanach", for- 
manl le message de navigation. 
[0061] Les couples de signaux LA et Q_A, LB et 

so q_b, LC et Q_C (ou leurs modules, ou le carre de ceux- 
ci), ainsi que le signal SIGN_IA_2K vont tous a Punite 
de gestion/decision 90 (liaisons non representees pour 
eviter de surcharger le dessin de la figure 2). En outre: 

55 - |es signaux LA et Q_A vont a un etage de boucle 
de porteuse 40, qui est de preference du type PLL 
(phase lock loop), et dont la sortie DVCa (ecart de 
vitesse porteuse), sur 16 bits, va, dans la phase de 
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poursuite du signal, d'une part vers un oscillateur 
de porteuse 71, a commande numerique, d'autre 
part vers un etage 59 qui definit un ecart de vitesse 
code DVCo, que Ton decrira plus loin ; 
les sorties LB et Q_B du demodulateur de porteuse 
32B de la voie B aboutissent a un etage de boucle 
de code 50, de preference du type DLL ("delay lock 
loop"). 

[0062] L'oscillateur de porteuse 71 opere a la caden- 
ce H100K.(100 kHz). Sa sortie P_PHI, representative 
de la phase de porteuse, sur 7 bits, alimente un gene- 
rates sinus-cosinus 75, egalement pulse a la cadence 
H100K, et fournissant sur 6 bits les signaux respectifs 
COSinus et SINus de la phase de porteuse, qui vont 
vers les demodulateurs de porteuse 32 deja cites. 
[0063] A partir du signal DVCa, l'oscillateur de porteu- 
se 71 ajuste sa phase P_PHI, mais aussi, chaque fois 
que celle-ci change modulo 2 n, en fournissant une in- 
formation d'increment ou decrement d'un tour de phase 
a un circuit 72 qui accumule alors, sur 16 bits, le nombre 
de tours P_DF que fait la porteuse, c'est a dire la partie 
entiere de sa phase totale. La partie fractionnaire de la 
phase totale, sur 16 bits, est la phase de porteuse 
P_PHI precitee. L'ensemble P_DF et P_PHI va a I'unite 
de gestion/decision 90. 

[0064] L'oscillateur de code 81, egalement un oscilla- 
teur a commande numerique, est pulse par une horloge 
a la frequence H25M (25 MHz). Sa sortie, sur 1 bit ali- 
mente un generateur de code 85. IL fournit Egalement 
une information de Iraction de chip C_PHI (ou "phase 
code"), sur 16 bits, qui est retourneea I'unit6de gestion- 
decision 90 (exploitation). 

[0065] La commande numerique de l'oscillateur de 
code 81 est definie sur 32 bits par la sortie d'un addi- 
tionneur 80, dont une entree recoit un signal CV, a une 
frequence de 1 ,023 MHz, qui est la cadence de base du 
code pseudo-aleatoire. L'autre entree de cet addition- 
neur 80 recoit une correction (aide en vitesse), prove- 
nant de I'unite 59 d6ja mentionnee. Cette correction 
DVCo correspond a la sortie de la boucle de code DLL 
50, "aidee" par la sortie de la boucle de porteuse PLL 
40. L'unite 59 permet de convertir la correction DVCa 
dans la meme unite de temps que celle de la correction 
DVCo. Le rappon de conversion K de I'unite 59 corres- 
pond au rapport des frequences respectives code et 
porteuse, soit 1540 (1575,45 divise par 1,023). 
[0066] Pour sa part, le generateur de code 85 recoil 
une information de calage temporel en provenance 
d'une table 84, rempliesous lecontrolede I'unite) de ges- 
tion-decision 90 en fonction du numero du code PRN 
relatif au satellite concern^. 

[0067] II lournit un code de base C/A, relatif a I'emet- 
teur concerne (frequ nee porteuse), sur une sortie al- 
lant vers I'unite de format de code 89. La partie fraction- 
naire C_PHI_FRAC de I'information "phase code" est 
lue en sortie de l'oscillateur de code, comme deja decrit. 
En outre, le generateur de code 85 fournit une informa- 



tion dite "registre G1 \ representative de son calage pro- 
pre, et qui aboutit a I'unite 86. Celle-ci convertit I'infor- 
mation du registre G1 en un N° ou rang de chip (1 parmi 
1023), qui represente la partie entiere C_PHI_INT de 

s I'information appelee "phase code" C_PHI. 

[0068] A partir du code de base C/A genere par le cir- 
cuit 85, on va fabriquer les repliques A, B : C de ce code 
correspondant aux trois voies de correlation decrites 
precedemment. Ces repliques permettent de piloter les 

io correlateurs correspondents 31 1 A a 31 1 C; Tune, voie A, 
est dite ponctuelle; les deux autres B et C sont dites E- 
L ("Early minus Late"), avec des espacements diffe- 
rents. 

[0069] Cette generation des repliques est realisee 

is dans I'unite de formatage des codes 89. Le lonctionne- 
ment de ce formatage correspond a celui d'une ligne a 
retard dont le schema est donn6 en figure 5. 
[0070] Un registre a decalage 890, de 63 cases (- 
31.0.+31), est aliments par le code issu du g6n6rateur 

20 de code 85, 6chantillonne a 25 MHz. Chaque case du 
registre 890 represente done une p^riode de cette hor- 
loge a 25 MHz, soit 40 ns (nanosecondes), ou encore ^ 
1/25 chip. La sortie de 5 cases parmi les 63 possibles^ 
va permettre de generer les repliques A, B et C. 

25 [0071] Une sortie pre-programmee 891, issue de la- 
case centrale 0, representee Pentr6e de la voie ponc- 
tuelle centred, ou voie A, codee sur 1 bit. 
[0072] Les autres sorties 892 E et 892 L, ainsi que 
893 E et 893 L, sont programmabies deux a deux sy- 

30 metriquement par rapport a la case 0. Elles sont com- 
binees dans des operateurs logiques 894 et 895, qui ef- 
fectuent une soustraction logique et une division par 
deux, et fournissent respectivement les entrees des cor- 
relateurs des voies B et C, codees sur 2 bits. 

35 [0073] Ainsi, la replique locale de la voie B correspond 
a une position des sorties 892 E et 892 L a +/- x fractions 
de 1/25 chip, d'ou un ecartement ou espacement de 2.x/ 
25 chip (entre les deux codes Avance E et Retard L dont 
on fait la difference). 

40 [0074] De son c6te\ la replique locale de la voie C cor- 
respond a une position des sorties 893 E et 893 L a +/- 
y fractions de 1/25 chip, d'ou un ecartement ou espace- 
ment de 2 y/25 chip. 

[0075] Les valeurs x et y sont programmabies par les 
45 entrees x et y du bloc 89 de la figure 2 en provenance 
de I'unite de gestion-decision 90, avec: 

x variable de 1 a 15, cod6e sur 4 bits, et 
y variable de 1 a 32, codee sur 5 bits et y supe>ieur 
50 & x . 

[0076] La voie C peut etre formee soit par la sortie, 
dite C2, de 895 (de la forme 893 E moins 893 L), soit 
par la sortie de 893 E ou bien. celle de 893 L (sorties 
55 indifferemment notees C1). 

[0077] On notera que les codes A, B et C2 sont definis 

conjointement, et avec le meme centre. 

[0078] Comm indique plus haut, lorsque Ton correl 
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le signal re9u avec une replique locale du code PRN 
attendu (voie A), la reponse du correlateur possede une 
allure triangulaire, avec un maximum au synchronisme, 
et la largeur de la base du triangle est de +/- Dmc, ou 
+/- 1 chip. 

[0079] Lorsque Ton correle le signal recu avec la dif- 
ference entre deux repliques locales, espacees, du co- 
de PRN attendu (repliques "differentielles" Bou C2), le 
resultat de la correlation (elle aussi "differentielle") cor- 
respond a ce que donnerait la difference entre deux cor- 
relations faites respectivement avec les deux compo- 
santes de la replique diflerentielle: par exemple, la sortie 
du correlateur de la voie B avec la replique locale (892 
E - 892 L) correspond a la difference entre les sorties 
de deux correlations faites respectivement avec les re- 
pliques locales disponibles en 892 E et 892 L. L'ecart 
temporel entre 892 E et892 L (par exemple) est nomme 
espacement, ou fenetre d'ecartom6trie, ou encore ecar- 
tometre. Les deux correlations dont on construit ainsi la 
difference ont chacune une reponse triangulaire de lar- 
geur +/- 1 chip a la base. Dans la mesure ou I'espace- 
ment est inferieur a 2 chips, elles vont se composer. Ain- 
si, la figure 9A illustre la reponse de la voie B, avec un 
espacement 2x = 0,2 chip (cette valeur de x est prise a 
titre d'exemple), tandis que la figure 9B illustre la repon- 
se de la voie C, avec un espacement 2y = 1 chip. Sur 
ces figures, A symbolise la duree d'un chip, notee Dmc 
(=2.v) dans la description, alors que le coefficient A est 
une grandeur liee a Pamplitude du signal recu. Dans les 
deux cas, cette composition fournit, au centre, une zone 
de reponse sensiblement lineaire, qui passe par zero 
au point milieu de la fenetre de correlation diflerentielle. 
C'est cette zone qui permet le pilotage d'un asservisse- 
ment temporel, comme le fait la voie B. Normalement 
(Figure 9A), on observe de part et d'autre de cette zone 
centrale deux paliers, de duree - d + Dmc, ou d est I'es- 
pacement de I'ecartometre ; la duree est done 2(v-x) sur 
la figure 9A. Sur la figure 9B, le palier se reduit a un 
point car 2(v-y) = 0, du fait que v=y. 
[0080] A I'initialisation du recepteur, une phase con- 
nue dite acquisition ou accrochage, permet d'amener le 
code ponctuel de la voie A suffisamment pres du syn- 
chronisme avec le signal recu pour obtenir une reponse 
(un pic) dans les correlateurs 311 A et 31 2A, et, par con- 
sequent, des signaux l_A et Q_A se trouvant au dessus 
du bruit. Ceci peut se faire par recherche systematique, 
et/ou en tenant compte d'informations qui sont deja dis- 
ponibles, de source externe ou bien du fail du fonction- 
nement anterieur du recepteur GPS (re-accrochage). 
On peut diminuer le temps de recherche systematique 
en utilisant en parallele le code A, ainsi que I'une et/ou 
Tautre des versions 893 E et 893 L du code C (eventuel- 
lement amenager la voie B, pour lui faire jouer a ce stade 
le meme role que la voie C). Cett r cherche syst 'ma- 
tique p ut se faire de toutes les manieres connues de 
rhomme du metier. 

[0081] Ensuite, on entre en phase de poursuite, avec 
les codes A, B et C2. Un signal D_B. dont on reparlera 



plus loin en reference a la figure 6 (sortie de I'etage 50 
d'ecartomet'ie avec boucle a verrouillage de phase 
PLL), pilots ia boucle de code jusqu'a obtenir un calage 
tempore! a +/- un dixieme de chip, ou mieux, par la voie 
5 B. Le positionnement grossier de la boucle de code (a 
+/- 1 chip) peut etre trouve a I'aide d'un signal d'ecart 
D_C, dont on reparlera egalement plus loin en reference 
a la figure 6 (sortie de Petage 60 d'ecartometrie). 
[0082] En regime permanent, I'asservissement de 
10 poursuite s'effectue avec les codes A et B. Sous I'effet 
de I'asservissemen1 de code de la voie B, la reponse de 
la voie A (sortie des accumulateurs 321 A et 331 A) est 
a un niveau d'autant plus eleve que le synchronisme en- 
tre le code recu et le code A, ou code ponctuel, est plus 
is precis. Ensuite, les multipiieurs 322A, 323A, d'une part, 
ainsi que 332A, 333A, d'autre part realisent la demodu- 
lation par la phase de ia porteuse. Les recombinaisons 
realisees par les organes 325 et 335 donnent l'ecart de 
phase, en ses composantes exactes, compte-tenu de 
20 ia separation en deux voies complexes I el Q. Apres ac- 
cumulation en 327A et 337A, les signaux LA et Q_A 
attaquent I'etage PLL de boucle de porteuse 40, qui est 
par exemple une boucle du 2° ordre, avec une bande 
passante de 30 Hz. Sa sortie est representative de 
25 l'ecart de phase de porteuse, et sert de commande nu- 
merique pour ajuster Poscillateur 71, dont la frequence 
de base est H100k (100 kHz). Cet oscillateur 71 forme 
memoire de phase totale pour la porteuse, avec: 

30 - une sortie P_PHl de phase modulo 2 pi, appliquee 
au generateur de cosinus et sinus 75 (sur 7 bits), et 
a I'unite de gestion/decision 90 (sur 16 bits, pour 
preserver la finesse du decalage Doppler), et 
une sortie P_DF de nombre de tours, appliquee a 

35 I'unite de gestion/decision 90. Ce nombre de tours 
represente une variation de frequence. 

[0083] Rapportees a la frequence de 100 kHz : ces 
deux grandeurs P_DF et P_PHI definissent ensemble 

40 une mesure fine de l'ecart temporel sur la voie A, en 
utilisant la frequence porteuse, au voisinage du syn- 
chronisme sur le code recu. Cette mesure fine est ap- 
pelee pseudo-distance porteuse (ou "pseudo range" 
porteuse). En phase d'accrochage (acquisition) en bou- 

45 cle ouverte, la valeur de consigne de I'oscillateu r de por- 
teuse est fournie par I'unite de gestion/decision 90 en 
fonction d'un Doppler estime ("Dop. Port"). 
[0084] En phase de poursuite en boucle fermee, la 
grandeur d'ecart correspondante en sortie de I'etage 

so PLL 40 est appliquee a I'oscillateur de porteuse, ainsi 
qu'au registre 59 eta I'additionneur numerique 80, pour 
"aider en Vitesse" la boucle de code asservie sur les 
composantes en phase LB et en quadrature Q_B de la 
voie B. La boucle de code st par exemple une boucle 

55 du 2° ordre avec une bande passante de 1 Hz, aidee en 
Vitesse par une sortie de ia boucle de poursuite porteu- 
se. 

[0085] En bref, dans un etat de poursuite stabilisee 
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[0086] Dans ce qui precede, la phase estimee de la 
porteuse <J>j est une phase totale, c'est-a-dire qu'elle 
peut comporter un decalage de frequence d'un certain 
nombre de tours, comrne on I'a vu. II est done plus sim- 
ple de la decomposer en une frequence f j, et une phase 
proprement dite, c'est-a-dire modulo 2n, notee (fc. 
[0087] Ces parametres de retard estime de transmis- 
sion du code entre le satellite et le recepteur et de fre- 
quence Doppler (avec phase) estimee dont la porteuse 
du signal est affectee, doivent permettre de piloier les 
correlations afin : 

de detecter la presence d'un signal, ce que fait la 
voie A A, afin de valider les estimations grossieres de 
t et f utilisees dans la phase de recherche, 

toujours a I'aide de la voie A, de commander la bou- 
cle de poursuite porteuse, qui est d'abord accro- 
chee en frequence, puis en phase, des que Pecart 
de frequence entre le decalage Doppler reel et f de- 
vient suffisamment faible : 

simultanement a I'accrochage de la boucle de pour- 
suite porteuse, ia boucle de poursuite code de la 
voie B est enclenchee; cette boucle permet de ren- 
dre plus fine Pestimation t de la distance satellite/ 
recepteur. Enfin, cette estimation plus fine x sera a 
son tour lissee par I'analyse de phase faite a I'aide 
de la boucle de porteuse. On entend ici par lisser, 
par exemple le fait que Ton utilise les valeurs de 
phase recalees par les valeurs moyennes de dis- 
tances. 

[0088] La boucle de poursuite porteuse elabore diffe- 
rents signaux de service relatifs a son etat, et transmis 
a Punite de gestion/decision 90, en particulier: 

un signal logique IPPP (pour Indice de Poursuite 
Porteuse en Phase) indiquant que la boucle de por- 
teuse est passee du mode de poursuite en frequen- 
ce au mode de poursuite en phase, (au besoin 
qu'elle reste de facon stable en poursuite en pha- 
se), 

un signal logique IRBP (pour Indice de Rejection 
par la Boucle de Porteuse) indiquant que la rejec- 
tion effectuee par la boucle de porteuse est corr c- 
te. Cet indice est fonde, de preference, sur I'analy- 
se, selon des criter s predetermines, des niveaux 
dans les differentes cases disponibles en sortie de 
la FFT, au niveau du discriminateur de la boucle de 



poursuite porteuse. Le critere le plus simple,, non 
limitatif, est que Penergie soit essentiellemerit si- 
tuee dans la case centrale (ou generalement la ca- 
se principale), et non dans les cases adjacentes a 
celle-ci; 

un signal logique IQPP (pour Indice de Qualite de 
Phase Porteuse). Cet indice est par exemple fonde 
sur le fait que le bruit du comparateur de phase de 
la boucle 40 est ou non acceptable, 
eventuellement, un signal logique IQPC (pour Indi- 
ce de Qualite de Poursuite Code), indiquant que la 
boucle de code ne trouve pas son equilibre (situa- 
tion en principe impossible lorsque la voie B alimen- 
te la poursuite distance - sur le code -, via le produit 
scalaire et la normalisation obtenus a partir de la 
voie A, comme e'est ici le cas), 
eventuellement aussi, un signal logique ITPP (pour 
Indice de Trainee en Poursuite Porteuse), indiquant 
que la boucle de porteuse peine a suivre une forte 
variation de frequence du signal incident, due en 
principe a une forte variation de la composante 
d'accel^ration le long de la ligne de visee. 

[0089] Les recepleurs GPS connus comportent deja 
25 des moyens de decision pour ecarter une ligne de visee 
qui fournit une mesure douteuse. La decision d'invalider. 
une ligne de visee se fonde par exemple sur le fait que 
la boucle de porteuse ne demeure pas d'une maniere 
stable en mode de poursuite de phase (signal logique 
so ippp), ou bien que la boucle de code ne trouve pas son 
equilibre (signal IQPC). De telles dispositions existent 
dans beaucoup de recepteurs GPS, pour invalider des 
lignes de visee et/ou pour indiquer a Pusager qu'une li- 
gne de visee est douteuse. Le choix des signaux utilises 
35 comme base pour cela depend largement de la concep- 
tion des boucles de poursuite de porteuse et de code. 
Si la ligne de visee est invalidee, le recepteur ne valide 
pas les estimees correspondantes, et travaille seule- 
ment avec les estimees disponibles pour les autres li- 
40 gnes de vis§e. II peut tenter une nouvelle phase d'ac- 
crochage ou re-accrochage avec le meme satellite, ou 
plutot rechercher un autre satellite plus propice. 
[0090] II apparaTt done que le recepteur est sensible : 

45 a) au niveau moyen d'^nergie du signal el a la loi 
de fluctuation d'amplitude pour la detection, 

b) a ia repartition spectrale de Penergie pour ia pour- 
suite porteuse, 

c) a la repartition tempore He de Penergie pour la 
50 poursuite code. 

[0091] La Demanderesse a observe que la presence 
de multiples trajets indesirables agit sur ces trois carac- 
teristiqu s et peut fauss r I s mesures generees par I 
55 recepteur pour un satellite donne. on sait que, pour des 
signaux sans interference, t pour un GPS absolu, I'er- 
reur de position peut aller jusqu'a 175 metres (en tri- 
dimensionnel). Cette erreur tombe de 1 a 10 metres en 



(en Pabsence de multiples trajets): 

la boucle de porteuse est accrochee en frequence 
et en phase, et 

la boucle de code est egalement accrochee, avec 
un signal d'ecart en principe proche de zero en sor- 
tie de la voie B. 
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GPS differential ; en interferometrie differentielle (geo- 
desie, determination d'attitude), I'erreurdevientde quel- 
ques centimetres, toujours en mode tri-dimensionnel. 
Dans ces applications, ou la classe de precision recher- 
chee est nettement meilleure, lorsqu'on utilise par 
exemple un GPS differentiel en code ou en phase, ou 
bien en interferometrie a base courte, le phenomene de 
multiples trajets devient meme un obstacle majeur, qui 
compromet la precision attendue du recepteur GPS. 
[0092] La Demanderesse a egalement remarque que 
la probability d'apparition simultanee de trajets multi- 
ples sur differents canaux est faible. En outre, du fait 
que Ton travaille de lacon generale en redondance des 
mesures disponibles, avec plusieurs satellites, il est ge- 
neralement possible de continuer le calcul de la solution 
de navigation en eliminant une, voire deux mesures 
douteuses. 

[0093] Dans cet esprit, la Demanderesse a recherche 
des techniques qui permettent d'identifier ou detecter 
les mesures douteuses en raison de multiples trajets, 
sans chercher a priori a en deduire des corrections de 
ces mesures douteuses. 

[0094] Les techniques ici proposees font intervenir 
une double analyse, a savoir une analyse frequentielle, 
effectuee au niveau de la boucle de porteuse, et une 
analyse temporelle, effectuee en liaison avec la boucle 
de code, et a I'aide d'une voie de correlation supplemen- 
tal. 

[0095] En bref, et comme on le verra, I'analyse fre- 
quentielle vise d'abord a detecter ies trajets multiples 
dont la composante Doppler se demarque de celle du 
trajet direct, alors que I'analyse temporelle vient la com- 
pleter en traitant la detection des trajets multiples dont 
I'ecart Doppler avec le trajet direct est faible ou tres fai- 
ble. 

[0096] On decrira tout d'abord I'analyse frequentielle. 
Cette technique offre une synergie particulierement in- 
teressante avec ('analyse temporelle, mais conserve 
aussi un interet propre, independamment de cette ana- 
lyse temporelle. 

[0097] Sur la figure 6, le bloc 38A designe un circuit 
de transformation de Fourier rapide, ou FFT, avec pour 
signaux d'entree l_A et Q_A. 

[0098] Cette translormation de Fourier rapide s'ettec- 
tue sur 512 points a 1 milliseconde, par exemple. Le cir- 
cuit d'analyse 39A recevant la sortie de la transforma- 
tion de Fourier va effectuer differentes operations que 
Ton delaillera plus loin, en liaison avec I'unile de gestion/ 
decision 90 (en pratique, le circuit 39A est incorpore a 
I'unite 90). 

[0099] Le fonctionnement du circuit d'analyse 39A est 
conforme au mecanisme de surveillance frequentielle 
illustre sur la figure 7. Apres le bloc d'entree 2000, la 
surveillance frequentielle comprend un premier test 
2001, qui consiste a examiner si le niveau de signal P 
a I'entree est normal, c'est-a-dire egal au bruit thermique 
KTBF, a une tolerance pres. Si cette condition n'est pas 
verifiee, il s'agit d'une anomalie generale - un dysfonc- 



tionnement serieux susceptible de porter sur differentes 
parties du recepteur - done en principe non diagnosti- 
cable en termes de multiples trajets. 
[0100] La suite du mecanisme fait intervenir des si- 

5 gnaux logiques, tels IQPR emanant notamment de la 
boucle de poursuite porteuse 40. 
[0101] Si I'indice IQPP est faux (etat 2010), indiquant 
par exemple que le bruit du comparateur de phase est 
excessif, alors l'etape2012 invalide la Vitesse, ladistan- 

io ce et la phase, pour le canal considere. 

[0102] Cette situation correspond au cas de la figure 
8A (dans ce qui suit on se referera sans distinction aux 
figures 8i-1 et 8i-2, ou \- A & G), ou Ton trouve en meme 
temps dans la case centrale de la sortie de la FFT (voir 

is plus particulierement les figures 8i-2) la contribution 
d'un signal direct, ainsi que d'un signal de reflexion spe- 
culate fort. La °ligne de visee" ne peut pas etre definie 
d'une maniere stable. S'agissant d'un recepteur GPS 
installe a bord d'un aeronef , cette situation peut se pre- 

20 senter notamment lors d'un virage au-dessus d'un lac, 
ou d'une mer d'huile, a basse ou moyenne altitude, ce 
qui se traduit par une reflexion speculate proche du si- 
gnal direct. II peut eventuellement s'y ajouter un,fond 
spectral ou des raies fortes supplementaires. 

25 [0103] La suite du mecanisme de la figure 7 fait inter- 
venir: 

les differentes cases de la FFT, designees par un 
indice i, egal a c pour la case centrale, ou se trouve 

30 le signal direct correspondant a la ligne de visee, 
un domaine dit "bande de securite", centre sur la 
case centrale, et dont ia largeur de ce domaine de. 
securite correspond avantageusement a une rejec- 
tion significative par la boucle de poursuite de por^ 

35 teuse, typiquement -10 dB; on note PBSj la puis- 
sance recue dans la case i de la bande de securite, 
avec i compris entre c-s et c+s, en excluant en prin- 
cipe sa case centrale (i * c); 
un domaine dit "bandes laterales" correspondant au 

40 reste des cases de la FFT, soit j < c-s ou j > c+s; 
PBLj indique maintenant la puissance (ou niveau) 
obtenu dans la case j de I'une des deux bandes la- 
terales, situees de part et d'autre de cette bande de 
securite, 

45 - des grandeurs de seuils BT0 et BT1, definies par 
rapport au bruit thermique, avec BT1 > BT0. Typi- 
quement, BT0 est defini en fonction du nombre de 
points de la FFT, et eventuellement du fait que Ton 
aura fait la moyenne de plusieurs carres de modu- 

50 ies de FFT successives. BT1 est choisi en vue d'une 
definition nette d'un "fort niveau" par rapport a BT0, 
comme on le verra plus loin. Dans certains cas, il 
peut etre avantageux que ces grandeurs de seuil 
BT0 t BT1 soient r ndues dynamiques. 

55 

[01 04] Le test 2020 (figure 7) porte, pour chacune des 
cases exceptee la case centrale d'indice c, sur le fait 
que la puissance PBSj franchit le seuil BT0 dans I'une 
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au moins des cases FFT de la bande de securit ' (voir 
par exemple figures 8B et 8D). Si aucun franchissement 
de ce genre n'existe (en dehors de la case centrale), 
I'etape 2021 constate que les donnees de phase et de 
vitesse sont valides, !e statut de la donn§e de distance 
restant a determiner. Si, au contraire, PBSj depasse le 
seuil BTO dans au moins Tune des cases consid6r6es 
de la FFT (hors la case centrale), I'etape 2022 etablit 
que les donnees de phase sont invalides et que les don- 
nees de distance sont en principe invalides, le sort des 
donnees de vitesses restant a preciser. 
[0105] On considere maintenant la sortie de I'etape 

2021. Le,test 2031 recherche aussi si, dans au moins 
l'une des cases de bandes laterales, PBLj depasse le 
seuil BTO. Si la reponse est non, I'etape 2032 etablit que 
les donnees de distance sont valides, et I'on passe di- 
rectement a I'etape de fin 2090. En effet, dans ce cas, 
toutes les cases de la FFT sauf la case centrale sont en 
dessous du seuil BTO. 

[0106] Si le test 2031 indique que PBLj depasse le 
niveau BTO dans Tune au moins des cases des bandes 
laterales (figure 8C), alors I'etape 2033 etablit que les 
donnees de distance sont en principe invalides, sous 
reserve de Teyentuelle correction dont on parlera ci- 
apres. 

[01 07] Cela peut corresponds a une reflexion specu- 
late forte isolee dans une case FFT des bandes latera- 
les (Fig. 8C). Ce cas se rencontre lorsqu'un aeronef ef- 
fectue un virage au dessus du sol, lequel comporte un 
objet fixe reffechissant, et produisant un defilement 
Doppler, vu de l'a6ronef. 

[0108] Si, en revanche, ce defilement Doppler est fai- 
ble, la raie de reflexion speculate forte se trouve alors 
dans la bande de securite (Fig. 8D), et a engendre un 
passage a l'6tape 2022. 

[0109] On retourne maintenant a la sortie de I'etape 

2022. Un test 2041 compare le niveau dans chacune 
des cases de la bande de securite, exceptee la case 
centrale, au seuil BT1 qui est superieur a BTO. Ce test 
recherche done s'il y a un franchissement de fort niveau 
du seuil dans la bande de securite hormis la case cen- 
trale. 

[0110] Si la reponse est oui, I'etape 2042 determine 
que les donnees de vitesse sont egalement invalides, 
et I'on passe a I'etape de fin 2090. 
[0111] Si la reponse est non, I'etape 2043 admet que 
les donn6es de vitesse sont valides : dans la mesure ou 
le signal logique IRBP indique que la rejection effectuee 
par la boucle de porteuse est correcte. 
[0112] On considere le cas ou un emetteur perturba- 
teur, proche du recepteur, provoque une interference, 
mais sans que la puissance a ('entree du recepteur ne 
s'en trouve sensiblement modified. Si cette emission 
parasite est laterale (Fig. 8E), lie depasse le seuil BTO, 
mais non I seuil BT1 , dans la bande de securite, et n'est 
pas assimilable a une reflexion diffuse; d'ou le passage 
par I'etape 2043 dans ce cas. 

[0113] Si au contraire remission parasite depasse le 



seuil BT1 dans la bande de securite (Fig. 8F), alors on 
passe par I'etape 2042. 

[0114] Apartirdes etapes 2033 ou 2043, on peut pas- 
ser directement a I'etape de fin 2090. La figure 7 illustre 
5 une variante optionnelle, sur laquelle on reviendra plus 
loin. 

[0115] On considere maintenant ranalysetemporelle, 
en mode poursuite. 

[0116] Sur la figure 6, la sortie de I'etage 50 d'ecarto- 
io metrie avec boucle a verrouillage de phase (PLL) de la 
boucle de code foumit une grandeur d'Scart D_B. 
[0117] La voie C recoit la replique C2 (893 E moins 
893 L); elle a done le meme centre que la voie B, mais 
un espacement plus grand, I'ecart temporel entre 893 E 
is et 893 L etant superieur a celui entre 892 E et 892 L 
Une unite d'ecartometrie 60 (identique a Punite 50, a 
I'exception de I'espacement entre les portes qui subdi- 
visent I'espacement) construit une grandeur d'ecart 
D_C, a partir des composantes en phase l_C et en qua- 
20 drature Q_C disponibles en sortie de la voie C. 

[0118] En outre, un soustracteur 61 construit la diffe- 
rence 

25 Z = D_B - D_C. 

[01 1 9] L'analyse temporelle comprend la surveillance 
du module du signal D_C (ou de son carr6), et/ou la sur- 
veillance du module du signal Z (ou de son carre). La 

30 figure 10 montre que ces modules iD_C! et !zi sont com- 
pares (901,902) a des seuiis respectifs Mc et M z . Un 
module de decision 910 utilise le depassement des 
seuiis pour construire une indication de la presence de 
multiples trajets, comme on le verra ci-apres. Bien que 

3S Ton prefere actuellement utiliser les deux signaux D_C 
et Z, il est envisageable d'utiliser un seul de ces deux 
signaux, du moins pour certaines applications. En pra- 
tique, les elements 901, 902 et 910 font partie de I'unite 
de gestion/decision 90. 

40 [0120] En outre, I'invention n'exclut pas la surveillan- 
ce du module du signal D_B (ou de son carre), en par- 
ticulier dans un recepteur qui n'aurait pas de voie C. La 
largeur d'ecartometrie de la voie B peut alors etre aug- 
mentee. 

45 [0121] On rappelle que les signaux D_B et D_C re- 
sultent de I'integration sur un nombre important de cor- 
relations individuelles. Cette integration peut aller de 
quelques centaines de millisecondes a quelques secon- 
des, suivant la dynamique du phenomene de trajets 

so multiples pour I'application consideree. De plus, I'espa- 
cement de la voie C etant nettement plus large que celui 
de la voie B, ia contribution du bruit dans le signal D_C 
est nettement superieure a celle du bruit dans le signal 
D_B. 

55 [0122] A titre non limitatif, on considere ci-apres que 
x et y sont choisis pour correspondre a des espace- 
ments de 0,2 chip et 1 chip, respectivement, pour les 
voies B et C. La valeur 2x peut descendre jusqu'a 0,1 
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chip environ; quant a y, on observe une attenuation as- 
sez rapide au dela de 1 chip. De preference, y est pris 
aussi grand que possible par rapport a x. 
[0123] La voie B alimente (a poursuite distance (sur 
le code), via le produit scalaire el la normalisation obte- 
nus a partir de la voie A. Apres un regime transitoire plus 
ou moins long, Passervissement r6alis6 par la boucle de 
code trouve en principe toujours une position d'equili- 
bre, du moins lorsqu'on travaille sur la base du produit 
scalaire. 

[01 24] Des lors que la fenetre ou ecartometre tempo- 
rel de la voie B, d'espacement 2.x etroit, recoit seule- 
ment le trajet direct, la grandeur D_B est a zero, abs- 
traction taite du bruit, lorsque Passervissement realist 
par la boucle de code est stabilise (comme rappele a 
I'interieur du bloc 50). Dans ces conditions, comme la 
voie C a le m§me centre, la grandeur D_C est ^gale- 
ment a z6ro, abstraction taite du bruit, mais sous reser- 
ve que la tenetre ou Ecartometre temporel de la voie C, 
d'espacement 2.y grand par rapport a 2.x, recoive seu- 
lement le trajet direct, elle aussi. 

[0125] La poursuite est dite canonique lorsqu'il n'y a 
qu'un trajet direct, et rien d'autre (au bruit pres) entre 
I'eYnetteur (le satellite) et le recepteur. Dans ce cas, le 
signal direct est associe a un retard t d (par rapport a la 
base de temps locale) qui est normalement voisin du 
centre x de la voie B. Ce signal aura les memes contri- 
butions sur la voie B et sur la voie C, car \\ entre a la fois 
dans les fenetres de ces deux voies. 
[01 26] Par contre, dans une situation de multiples tra- 
jets, le signal direct S D est accompagne d'un signal re- 
flechi S Rl au moins. On note t r le retard du trajet refle- 
chi, et d RD la difference x R - t D . Les travaux de la De- 
manderesse lui ont permis de faire plusieurs constata- 
tions: 

a) si deux signaux GPS portant le m§me code sa- 
tellite, par exemple Tun direct, Pautre reflechi, sont 
simultanement presents a I'interieur d'un meme 
ecartometre (voie B et/ou C) : celui-ci prend une po- 
sition d'equilibre x qui est differente a la fois du re- 
tard t d du trajet direct et du retard t r du trajet refle- 
chi. Le module de la sortie decet ecartometre, hors 
bruit, est d'autant plus different de zero que I'ecart 
de phase entre les signaux direct S D et reflechi S R 
est plus different de kn (k entier, et phases definies 
par rapport a la frequence porteuse, Doppler com- 
pris). Une premiere possibility de detection d'ano- 
malie (indication de la presence de multiples trajets) 
peut done etre fournie par le fait que le module de 
la sortie de cet Ecartometre depasse un seuil. Ce 
seuil tient compte du niveau de bruit attendu ou me- 
surd; il est sensiblement proportionnel au temps 
d'integration; il tient compte aussi d'un risque d'er- 
reur admis ("probabilite de fausse alarme" ou Pfa). 
A ce stade, les signaux D_B et/ou D_C sont utilisa- 
bles. 



b) si les deux signaux GPS sont a krc pres (s< 
phase ou en opposition de phase), la pc- 
d'equilibre x atteinte par Pecartometre est tellf , 

le module de sa sortie soit nul. Mais la valeur : : r 
s depend de I'ecart utilise entre les deux portes t :-t 
L de cet ecartometre. Si Ton utilise deux ecartome- 
tres, d'ecartement de portes differents, et que Ton 
pilote la poursuite avec Tecartometre etroit (comme 
e'est ici le cas), I'ecartometre large, de memecen- 
10 tre, fournit une sortie de module non nul, car sa pro- 
pre position d'equilibre i* ne correspond pas a x. 
Une seconde detection d'anomalie peut done se fai- 
re en utilisant deux ecartometres, Pun etroit pour la 
poursuite, Pautre large pour la surveillance. II a ega- 
15 lement ete observe que I'ecart entre les deux posi- 
tions d'equilibre t 1 et x est plus grand lorsque I'ecart 
de phase est n que lorsqu'il est nul (a 2kn pres). 
Ceci est particulierement interessant, puisque 
I'ecart n ou situation en opposition de phase est le 
20 cas ie plus genant. 

c) lorsque I'on n'est plus en mode stationnaire, par 
exemple au moment de I'apparition (rapide ou non) 
d'un signal reflechi S R , ou bien en presence d'un 

25 6cart de frequence significatif entre le signal direct 
S D et le signal refl6chi S R , il a et6 observe que le 
module de Z peut avantageusement remplacer ou 
completer Pindication donn6e par le module de 
D_C, du fait qu'il contient moins de bruit que celui- 
30 C i. 

[01 27] On precisera maintenant ce qui precede, en te- 
nant compte du bruit. 

[01 28] La Demanderesse a observe ceci: des lors que 
35 la boucle de code fonctionne de maniere normale et sta- 
bilisee sur le trajet direct seulement ("poursuite canoni- 
que' du code a Paide de la voie B), la contribution de 
Perreur residuelle de positionnement de la boucle de co- 
de dans !D_C! (le module du signal D_C) est nettement 
40 au dessous du bruit contenu dans D_C (typiquement 
dans un rapport de I'ordre de 10). Elle est done n6gli- 
geable. 

[0129] II est alors propose de fixer pourlD_Cl un seuil 
M c , dont le franchissement sera indicatif de la presence 

45 d'un trajet supplemental. Ce seuil peut §tre fixe par 
rapport au niveau de bruit attendu dans la voie C, et 
compte-tenu d'un risque d'erreur admis ('probabilite de 
fausse alarme" ou Pfa). Le seuil est proportionnel au 
temps d'integration. 

so [0130] Par exemple, pour 7 satellites disponibles, la 
Pfa est de 14%, si I'on admet de rejeter en moyenne 
une ligne de vis6e pour chaque jeu de 7 mesures simul- 
tanees utilisables pour calculer un positionnement. On 
peut alors en deduire par caicul numeriqu le seuil a 

55 appliquer a la voie C, compte tenu des parametres de 
fonctionnement du r'cepteur. 

[0131] Si !D_C! depasse M c , I'unite de gestion/deci- 
sion 90 va done decider d'invalider la ligne de visee en 
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cours. On peut faire de meme avec !D_Bl, pour la voie 
B, comme deja indiqu6. 

[0132] On verra maintenant que le signal Z contient 
d'interessantes informations sur I'existence et I'amplitu- 
de des signaux issus de trajets autres que le trajet direct. 
[0133] Dans le cas d'une poursuite canonique, la dif- 
ference Z ne comporte que du bruit. Mais la Demande- 
resse a observe que le bruit dans !z! = lDJ3 - D_Cl est 
inferieur au bruit dans !D_Cl, et ce typiquement d'envi- 
ron 20 % pour un rapport 5 entre y et x (y=5x). 
[0134] II est apparu que, dans la difference DJ3 - 
D_C, le bruit de D_B se soustrait partiellement du bruit 
de D_C: dans la partie commune des fenetres de la voie 
C et de la voie B, il existe une composante de bruit qui 
depend du positionnement et de Tecart, et est done dil- 
ferente dans Tune et I'autre voie, mais aussi un fonds 
tempore! de bruit, qui est le meme, et disparait par sous- 
traction. En d'autres termes, e'est leur espacement dif- 
ferent qui d£correle plus ou moins les sorties des deux 
correiateurs des voies B et C. 

[0135] Sur la figure 9C, I'abscisse represente I'ecart 
temporel d RD (sans preoccupation d'echelle), tandis 
que I'ordonnee represente !z!, a une echelle definie de 
la maniere suivante: 

AX = A R (1 - x/v) 

ou A R = N a r 0,6/o , a r etant I'amplitude reelle du signal 
reflechi S R , tel que recu, N le nombre de points d'echan- 
tillonnage (le produit de la frequence d'echantillonnage 
F e par le temps d'integration dt). et o I'ecart-type du bruit 
dans la bande utile en sortie de I'etage d'entree du re- 
cepteur. 

[0136] Sur la figure 9C, on considere en outre que le 
signal direct est toujours compris dans Tecartometre de 
la voie B, tandis que le signal reflechi va "balayer" I'in- 
tervalle de temps de I'ecartometre de la voie C. 
[0137] Si d RD est inferieur ou egal a +/- x, le signal 
reflechi (S R ) et le signal direct (S D ) auront les memes 
contributions dans la voie B et dans la voie C. La diffe- 
rence Z ne comporte done que du bruit. C'est ce qu'il- 
lustre la figure 9C entre les abscisses -x et +x. 
[0138] Pour le reste, la Demanderesse a observe 
que: 

lorsque d RD est compris entre>x et 2.v-y, il demeure 
une contribution croissante due au signal reflechi 
S R dans la voie B, laquelle, avec la contribution de 
S R dans la voie C, produit une croissance jusqu'a 
une valeur maximale +AX; 

lorsque d RD est compris entre 2.v-y et 2.v-x, la si- 
tuation est la meme que ci-dessus, sauf que la r6- 
sultante devient decroissante, desc ndant de +AX 
a +AX/2; 

lorsque d RD est compris entre 2.v-x et 2.v+x, il n'y 
a plus de contribution croissante du au signal re- 
flechi S R dans la voie B. Sa contribution dans la voie 



C demeure constante a +AX/2; 
lorsque d RD est compris entre 2.v+x et 2.v+y, la con- 
tribution due au signal reflechi S R dans la voie C 
decroit vers zero; 
5 . enfin, au dela de 2.v+y, le signal reflechi S R n'est 
plus vu par la voie C. 

[0139] La courbe est naturellement symetrique par 
rapport au point (0,0), pour les valeurs negatives de d RD 
10 (dans certaines applications, il peut arriver qu'un signal 
"refl§chi* - ou assimile, en principe un signal de leurrage 
- arrive avant le signal "direct"). 
[0140] Pour bien comprendre la figure 9C, il est ne- 
cessaire de tenir compte des elements suivants: 

15 

on a trace le module de la difference D_B - D_C, 
la poursuite distance sur la voie B s'effectue sur ce 
que "voit" la voie B de la resultante des signaux S D 
et S R , 

20 - on a suppose x D compris entre -x et + x, 

en outre, la representation de la figure 9C suppose 
que 2.y vaut 1 chip, soit y = v ; s'il en est autrement, 
rinlervalle de surveillance possible est diminue. 

25 [0141] II resulte de ce qui precede que la presence 
d'un trajet reflechi en plus du trajet direct, et a plus de x 
unites temporelles de ce trajet direct, se traduit dans le 
signal Z par une composante continue, qui depasse le 
niveau de bruit que contient ce signal. 

30 [0142] L'application d'un seuil M 2 convenablement 
choisi au signal Z permet d'identifier la presence de tra- 
jets multiples dans la voie C et/ou B, et d'en tirer les 
consequences sur le bon fonctionnement de la voie B, 
done de la boucle de code qu'elle pilote. Une facon - 

35 non limitative - de fixer ce seuil M z consiste a prendre 
le seuil M c prevu pour le signal D_C lui-meme, modifie 
a proportion du rapport des niveaux de bruit respectifs 
de D_C et de Z (80 %, si le bruit dans Z est 20 % sous 
le bruit dans D_C). Le signal !z! est done lui aussi un 

40 bon indicateur de ('existence de multiples trajets. 

[0143] Un autre facteur a prendre en compte est la 
"dynamique frequentielle" ou "dynamique porteuse" du 
signal incident. En cas de forte acceleration relative, il 
apparaTt un decalage de frequence de la porteuse. Les 

45 boucles de poursuite supportent cette dynamique fre- 
quentielle dans la limite de leur bande passante. II va 
en resulter cependant un retard ou "trainee" dans la r6- 
ponse de ces boucles. Et cette trainee cree une com- 
posante continue, qui apparait dans les signaux D_B et 

50 D_C, et peut faire par exemple que !D_C! depasse le 
seuil M c , en presence d'une forte acceleration relative 
lorsque les bandes passantes ne sont plus adaptees a 
la dynamique du recepteur, ce qui correspond a un mau- 
vais usage par I'utilisateur, et, alors qu'il n'y a pas de 

ss trajet multiple. On observera par contr que le signal Z 
est, de son cote, debarrasse de cette trainee, par diffe- 
rence. 

[0144] II est done toujours possible en considerant les 
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signaux D_C et Z, leurs seuils respectifs, et ies condi- 
tions de travail de la boucle de poursuite (en particulier 
le signal logique ITPP deja cite), d'obtenir une indication 
fiable de la presence de trajets multiples. On notera que 
le signal ITPP peut etre tourni par le fait que la detection 
frequentielle donne le signal principal non pas dans ia 
case centrale, mais dans I'une des cases voisines (ou 
a cheval sur des cases adjacentes voisines). La struc- 
ture a double ecartometre proposee est done particulie- 
rement interessante en elle-meme a cet effet. 
[0145] Les resultats de I'analyse frequentielle et/ou 
de I'analyse temporelle servent a ecarter selectivement 
les mesures (une "ligne de visee", ou plusieurs) identi- 
fies comme partiellement ou totaiement douteuses se- 
lon la presente invention. Les indicateurs fournis ici 
viennent done completer avantageusement les indica- 
teurs (independants des trajets multiples) qui sont -deja 
utilises dans les recepteurs connus pour ecarter une li- 
gne de visee fournissant une mesure douteuse. 
[0146] On considere maintenant le cas d'une invali- 
dation partielle ou par exemple seule la "vitesse" reste 
valide. Dans ce cas, il demeure possible d'utiliser des 
informations partiellement valides pour la resolution des 
equations de navigation, en fonction de la precision re- 
cherchee, sinon la ligne de visee est invalidee. 
[0147] Ajoutee a I'analyse frequentielle, I'analyse 
temporelle donne un moyen supplemental tres inte- 
ressant pour detecter la presence de trajets reflechis. 
Sur cette detection, on invalide toutes les mesures cor- 
respondent a la ligne de visee consideree. 
[0148] En effet, I'analyse frequentielle permet de se- 
parer des signaux reflechis du signal direct lorsque la 
distance du reflecteur a I'aerien du recepteur est varia- 
ble (ecart Doppler non nul). C'est par exemple le cas 
d'un avion en vol ascendant ou descendant par rapport 
a un sol reflechissant, ou bien le cas d'un vehicule ter- 
restre qui evolue a cote d'une surface potentiellement 
reflechissante, etnon paralleleau vecteur Vitesse du ve- 
hicule. 

[0149] En outre, I'analyse temporelle qui vient d'etre 
decrite donne une information la ou I'analyse frequen- 
tielle est peu efficace, c'est a dire a ecart Doppler faible 
ou nul entre le signal direct et le ou des signaux refle- 
chis. C'est par exemple le cas pour un aeronef en vol 
horizontal par rapport a un sol reflechissant, ou bien 
pour un recepteur stationnaire (fixe par rapport aux sur- 
faces reflechissantes). 

[01 50] On notera que les deux trailements sont entie- 
rement menes en parallele, et sans aucun multiplexage 
des correlations et autres asservissements. Us permet- 
tent en outre de traiter des multiples trajets aussi bien 
du type diff us que speculates, ce qu'aucun moyen con- 
nu n'a permis jusqu'a present. 

[0151] L'analyse fr 'quenti II est particulierement ef- 
ficace pour les multiples trajets diffus, et, a ce titre, in- 
teressante a elle seule. Plus precisement, les trajets dil- 
fus jouent sur la forme de la fonction de correlation des 
codes d'une maniere qui depend de la vitesse du porteur 



du recepteur, d'ou la necessite d'un temps d'analyse 
court. Dans ce cas. I'analyse frequentielle agil pleine- 
ment. 

[0152] Comparativement a I'analyse temporelle de- 

5 crite, il a deja ete propose d'analyser partiellement la 
forme de la courbe de correlation, notamment dans FR- 
A-2 698 966. Pour mener cette analyse plus profonde- 
ment, il faudrait soit prevoir un grand nombre de corre- 
lateurs, ce qui est lourd, soit, si Ton ne dispose que d'un 

10 nombre de correlateurs limite, decaler ceux-ci succes- 
sivement dans le temps, et prevoir un temps d'analyse 
eleve, en supposant en outre que la forme de la courbe 
de correlation est stable pendant tout ce temps d'ana- 
lyse, au moins. Comparativement, I'analyse par double 

is ecartometrie (voies B et C), telle que proposee, est sim- 
ple, immediate, et ne suppose pas une forme stable de 
la courbe de correlation sur un intervalle de temps long. 
[01 53] Comme on Pa vu, la combinaison de I'analyse 
frequentielle et de I'analyse temporelle permet de reje- 

20 ter les mesures de facon bien discriminatoire: 

s'il y a detection frequentielle de trajets multiples, 
speculates et/ou diffus, la mesure de distance sur 
la ligne de visee consideree est invalidee par pre- 

2S caution, 

si I'ecart Doppler est faible, on invalide les mesures 
de vitesse et phase issues de la boucle de porteuse, 
si I'ecart Doppler est grand, on valide par contre les 
mesures de vitesse et phase issues de ia boucle de 

30 porteuse, 

s'il y a detection de trajets multiples par I'analyse 
temporelle, on invalide toutes les mesures corres- 
pondantes, par precaution. 

3$ [0154] Bien que I'idee de base de I'invention soit de 
detecter la presence de trajets multiples, il demeure en- 
visageable de chercher a les corriger, du moins dans 
certains cas. Un exemple d'une telle correction option- 
nelle apparait sur la figure 7. 

40 [0155] Avec cette option de correction, a partir des 
etapes 2033 ou 2043, on aboutira au test 2051 , qui part 
d'une variable logique IPS (Indicateur de Possibilite de 
retour Sol), laquelle indique si I'altitude et I'elevation per- 
mettent ('observation des trajets multiples sur la ligne de 

45 visee en cours. Cette operation doit tenir compte de la 
relation suivante : 

2.H.sin0 < 5Dmin, 

so 

ou H est I'altitude, 0 est Tangle d'incidence des ondes 
sur le recepteur par rapport a I'horizontale ("elevation"), 
et 6Dmin represente un ecart plafond entre signal direct 
et signal refl 'chi, au-dela duquel la perturbation ne pro- 
55 voque plus de biais significatif en distance. Cette gran- 
deur SDmin est etablie en fonction de Pespacement du 
correlateur differentiel utilise dans la voie B; SDmin cor- 
respond typiquem nt a 250 metres dans I'exemple de- 
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crit. De leur cote, H et (-) sonl determinates, au besoin 
a I'aide des instruments cic bord d'un avion. Une valeur 
approchee suffisante de H peut etre obtenue en utilisant 
toutes les mesures (meme celles qui sont polluees par 
des trajets multiples). El (-) est calculable precisement, 
car le positionnement obtenu avec toutes les mesures 
sera faux au plus de 1 00 a 200 metres, ce qui est negli- 
geable angulairement. compte-tenu de la distance mi- 
nimale satellite-recepleur de 20 000 km. On notera in- 
cidemment que H est la hauteur par rapport au geoTde, 
et que Ton n'y accede qu'en ayant en memoire la hau- 
teur du geoide par rapport au WGS 84 dans la zone sur- 
volee, par exemple 50 metres dans I'ocean atlantique, 
en face des Landes. 

[0156] Si IPS est faux, la distance reste invalide com- 
me detini a I'etape 2033. et Ton passe a I'etape de tin 
2090. 

[0157] Dans le cas contraire. on passe au test suivant 
2053. Ce test porte sur I'existence de cases FFT avec 
f ranchissement du seuil BTO (abstraction faile de la ca- 
se centrale), ma is demeuranl isolees. Ce caractere iso- 
le est avantageusement determine par le fait que, pour 
chaque case franchie, les deux cases inferieures adja- 
centes et les deux cases superieures adjacentes ne font 
par contre pas I'objet d'un tranchisscmcnt du niveau de 
bruit BTO. Si la reponse est oui (case isolee), le signal 
I PR est vrai, et Ton passe a I'etape de fin 2090. Ceci 
correspond notamment au cas illustre en figure 8G, que 
Ton retrouve aussi en 8D et 8C. 

[0158] Dans le cas contraire, le test 2054 recherche 
s'il existe une seule serie de cases non isolees, faisant 
I'objet d'un f ranchissement du seuil BTO, et conservant 
un contraste important entre le signal direct et celui qui 
est contenu dans chacune de ces cases. En outre, la 
frequence centrale definissable sur I'ensemble de ces 
cases franchissant le seuil BTO doit demeurer relative- 
ment proche de la case centrale. Si tout cela est faux, 
on passe directement a I'etape de fin 2090, la distance 
restant invalide. 

[0159] Dans le cas contraire, il est envisageable a 
I'etape 2055 de corriger la distance d'apres I'integrale 
de I'energie dans cette serie ou "continuum" de cases 
entourant la case centrale, et possedant une frequence 
centrale relativement proche de celle de la case centra- 
le. Cette situation correspond au schema de la figure 
8B (le contraste est insuffisant sur la figure 8F). 
[0160] L'invention n'est pas limitee au mode de reali- 
sation decril, en particulier pour ce qui louche au recep- 
teur a trois voies a partir duquel l'invention est ici decrite. 
Un nombre superieur de voies n'est pas exclu. Dans cer- 
tains cas de figure, on pourrait envisager d'intervertir les 
roles des voies B et C. On pourrait aussi envisager que 
ces voies n'aient pas le mdme centre (entre elles, ou 
que la vol A), notamment pour tenir compte du fait que 
les traj ts reflechis naturels sont posterieurs au trajet 
direct. 

[0161] Defacon plusgenerale, l'invention peut s'eten- 
dre a d'autres systemes de radio-positionnement (ou 
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plus generalement de radio-navigation) que le GPS, 
qu'ils soient fondes sur des satellites, comme le syste- 
me GLONASS, ou bien terrestres. Dans le cas du sys- 
teme GLONASS, on sait que les satellites utilisent des 
s frequences porteuses differentes, modulees en bi-pha- 
se par un code pseudo-aleatoire qui leur est commun; 
I'homme du metier saura faire les adaptations requises 
pour la mise en oeuvre de la presente invention. 
[0162] De plus, on utilise ici les ondes pour localiser 
10 le recepteur; une variante d'application consiste a ree- 
mettre des ondes marquees temporellement en fonction 
des donnees de temps obtenues par la reception. On 
pourrait aussi appliquer l'invention par exemple a des 
stations de reference (dont la position est connue a prio- 
ns ri), permettant a des recepteurs GPS mobiles du voisi- 
nage de travailler en mode differentiel, bien que la con- 
trainte economique s'appiique beaucoup moins a la 
conception de telles stations de reference, et que 
d'autres moyens soient par consequent envisageables 
20 dans ce dernier cas. 



Revendications 

25 1. Dispositif recepteur de signaux radio-electriques 
comprenant une onde porteuse modulee munie 
d'un marqueur temporel a pseudo-bruit repetitif, du 
genre code pseudo-aleatoire, en particulier pour la 
radionavigation, 

30 le dispositif comportant une reception haute fre- 
quence, qui fournit un signal utile a plusieurs mo- 
dules a coincidence, dont chacun est affect e a un 
calage temporel sur un signal radio-electrique res- 
pectif choisi, et des moyens Sexploitation (90) des 

35 donnees de calage temporel fournies par ces mo- 
dules a coincidence, 

tandis que chaque module a coincidence com- 
prend: 

^0 - des moyens a memoire chronometrique 
(71-75,80-89), aptes a stocker, comme don- 
nees de calage temporel, une estimee de cala- 
ge de code (C_PHI) et une estimee d'etat de 
porteuse en frequence (P_DF) et phase 

4 5 (P_PHI), ainsi qu'a en tirer, avec une horloge 

locale, une image locale de porteuse et au 
moins une replique locale du code, pour le si- 
gnal radio-electrique concerne, et 
des asservissements de code et de porteuse 

50 (40,50,59), operant par correlation (311 A, 

312A;311B,312B) entre le signal utile et la re- 
plique locale de code, ainsi qu'a I'aide d'un si- 
gnal d'ecart de frequence et de phase (DVCa) 
entre la port use presente dans I signal util 

55 et ladit image de porteuse, afin de faire coin- 

cider la replique locale du cod avec le code 
recu, 
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caracterise en ce qu'il comprend des rnoyens de 
transformation de Fourier (38A) recevant le signal 
(DVCa) d'ecart de frequence et phase de porteuse, 
et des moyens de decision (39A,90) capables d'in- 
valider au moins partiellement les donnees de ca- 
lage temporel (C_PHI,P_DF,P_PHI), lorsque la 
transformation de Fourier indique la presence 
d'energie en dehors du voisinage d'une raie centra- 
le, correspondant a la coincidence, ce qui fournit 
une indication de presence de multiples trajets. 

2. Dispositif selon la revendication 1, caracterise en 
ce que les moyens de decision sont agences pour 
valider les estimees de frequence (P_DF) et de 
phase de porteuse (P_PHI), lorsque le niveau 
d'energie obtenu dans une bande de securite avoi- 
sinant la raie centrale, celle-ci etant exclue, teste 
inferieur (2020,2021) a un niveau de reference 
(BT0). 

3. Dispositif selon la revendication 2, caracterise en 
ce que les moyens de decision sont agences pour 
valider I'estim^ede calage decode (C_PH!) lorsque 
ie niveau d'energie obtenu dans des bandes latera- 
ls a la bande de securite reste inferieur 
(2031,2032) a un niveau de reference (BTO). 

4. Dispositif selon la revendication 2, caracterise en 
ce que les moyens de decision sont agences pour 
invalider au moins I'estimee de phase de porteuse 
(P__PHI), lorsque !e niveau d'energie obtenu dans 
ladite bande de securite, raie centrale exclue, de- 
passe (2020,2022) le niveau de reference (BTO). 

5. Dispositif selon la revendication 4, caracterise en 
ce que les moyens de decision sont agences pour 
invalider en outre I'estimee de frequence (P_DF) de 
porteuse et I'estimee de calage de code (C_PHI), 
lorsque le niveau d'energie obtenu dans la bande 
de securite depasse (2041 ,2042) un second niveau 
de reference (BT1), superieur au premier (BTO). 

6. Dispositif selon la revendication 4, caracterise en 
ce que les moyens de decision sont agences pour 
valider I'estimee de frequence (PJDF) de porteuse, 
lorsque le niveau d'energie obtenu dans la bande 
de securite ne depasse pas (2041 ,2043) le second 
niveau de reference (BT1), sous reserve du bon 
fonctionnement de I'asservissement de porteuse 
(IRBP). 

7. Dispositif selon I'une des revendications 2, 4 et 6, 
caracterise en ce que les moyens de decision sont 
agences pour valid r I'estimee d calage de code 
(C_PHI), lorsque les depassements dudit seuil 
(BTO) dans la bande de securite (2055) satisfont 
des conditions predetermines portant sur I'une au 
moins des proprtetes suivantes : continuity, centra- 



ge et non depassement du second seuil (BT1). 

8. Dispositif selon I'une des revendications preceden- 
tes, caracterise en ce que les moyens formant me- 
5 moire chronometrique comprennent un oscillateur 
de porteuse (71) pour definir I'image locale de por- 
teuse, un oscillateur de code (81 ), et un generateur 
(85-89) d'au moins une replique locale du code. 

10 9. Dispositif selon I'une des revendications preceden- 
tes, caracterise en ce que les asservissements de 
code et de porteuse (40,50,59) comprennent: 

une voie de correlation ponctuelle (A), compor- 
15 tant une demodulation selon ladite image de 

porteuse, et fournissant un signal de coinciden- 
ce (l_A,Q_A) entre la position temporelle du co- 
de dans le signal utile etson estimee de codage 
de code, 

20 - au moins une voie de correlation differenlielle 
(B,C), comportant une demodulation selon la- 
dite image de porteuse, et fournissant un signal 
d'ecart (D_B,D_C) entre la position temporelle 
du code dans le signal utile et son estimee de 
25 codage de code, 

le signal de coincidence (l_A,Q_A) servant de 
grandeur d'6cart pour une boucle de porteuse 
a verrouillage de phase, laquelle pilote le ge- 
nerateur (71 -75) d'image locale de porteuse, et 
30 - le signal d'ecart de code (D_B) servant de gran- 
deur d'ecart pour une boucle de code a ver- 
rouillage de retard, dont la sortie est combinee 
(59) a celle de la boucle de porteuse a ver- 
rouillage de phase, pour piloter le generateur 
35 (80-89) de repliques de code. 

10. Dispositif selon la revendication 9, caracterise en 
ce que les moyens de decision (90) sont agences 
pour etablir une autre indication de la presence de 

40 multiples trajets en utilisant le signal d'ecart (D_B, 
D_C) issu d'au moins une voie de correlation diffe- 
rentielte (B,C). 

11. Dispositif selon la revendication 10, caracterise en 
45 ce que le module a coincidence comprend au moins 

une seconde voie de correlation difference He (C), 
et en ce que les moyens de decision (90) sont agen- 
ces pour etablir encore une autre indication de la 
presence de multiples trajets en utilisant le signal 
so d'ecart (D_C) issu de cette seconde voie de corre- 
lation diff6rentielle (C). 

12. Dispositif selon la revendication 11, caracterise en 
ce que la seconde voie de correlation differ ntielle 

55 (C) possede le meme centre qu la premiere. 

13. Dispositif selon Tune des revendications 11 et 12, 
caracterise en ce que la seconde voie de correlation 
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differentielle possede un espacement (2.y) different 
de celui de la premiere (2.x). 

14. Dispositif selon la revendication 13, caracterise en 

ce que la premiere voie de correlation differentielle s 
(B) possede un espacement inferieur a 1 pas de co- 
de, et en ce que la seconde voie de correlation dif- 
ferentielle (C) possede un espacement superieur a 
celui de la premiere. 

10 

15. Dispositif selon la revendication 14, caracterise en 
ce que I'espacement (2.x) de la premiere voie de 
correlation differentielle (B) est de I'ordre du dixie- 
me d'un pas de code, tandis que I'espacement (2.y) 

de la seconde voie de correlation differentielle (C) * 5 
est de I'ordre d'un pas de code. 

16. Dispositif selon Tune des revendications 11 a 15, 
caracterise en ce que les espacements de la pre- 
miere et de la seconde voies de correlation diff6- 20 
rentielle sont programmables. 

17. Dispositif selon Tune des revendications 10 a 16, 
caracterise en ce que les moyens de decision (90) 
sont agences pour etablir ladite (ou lesdites) autre 2s 
(s) indication(s) de la presence de multiples trajets 
lorsque le signal d'ecart issu de la voie de correla- 
tion differentielle consideree depasse un seuil choi- 

si. 

30 

18. Dispositif selon Tune des revendications 11 a 16, 
caracterise en ce que les moyens de decision (90) 
sont agences pour etablir I'autre indication de la 
presence de multiples trajets en fonction d'une 
comparaison entre le signal d'ecart (D_C) issu de 35 
la seconde voie de correlation differentielle et le si- 
gnal d'ecart (D_B) issu de la premiere voie de cor- 
relation differentielle. 

19. Dispositif selon la revendication 18, caracterise en 40 
ce que les moyens de decision (90) sont agences 
pour etablir ladite autre indication de la presence 

de multiples trajets en fonction du fait que la diffe- 
rence (Z) entre le signal d'ecart issu de la seconde 
voie de correlation differentielle et le signal d'ecart 45 
issu de la premiere voie de correlation differentielle 
depasse un autre seuil choisL 

20. Dispositif selon I'une des revendications 18 et 19, 
caracterise en ce que ladite autre indication de la & 0 
presence de multiples trajets etablie en fonction de 

la difference (Z) tient compte d'un decentrage de la 
poursuite en frequence, indique par I'analyse fre- 
qu ntiell . 
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